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The Puzzles of RX J1856.5-3754: Neutron Star or Quark Star? 

J. E. Triimper a , V. Burwitz a , F. Haberl a , and V.E. Zavlin a 

a Max-Planck-Institut fiir extraterrestrische Physik, D-85741 Garching, Germany 

We discuss recent Chandra and XMM-Newton observations of the bright isolated neutron star RX J1856.5— 3754 
and suggest that the absence of any line features is due to effects of a high magnetic field strength (~ 10 13 G). 
Using different models for the temperature distribution across the neutron star surface assuming blackbody 
emission to fit the optical and X-ray spectrum and we derive a conservative lower limit of the "apparent" neutron 
star radius of 16.5km x (d/117pc). This corresponds to the radius for the "true" (de-redshifted) radius of 14km 
for a 1.4 Mq neutron star, indicating a stiff equation of state at high densities. A comparison of the result with 
mass-radius diagrams shows that quark stars and neutron stars with quark matter cores can be ruled out with 
high confidence. 



1. INTRODUCTION 

The excellent soft X-ray response of ROSAT 
pointed and all-sky observations pQ have led to 
the discovery of eight thermally emitting objects 
(listed in Table ^| having roughly the size of a 
neutron star but showing not any radio emis- 
sion |2 . 

Among these objects RX J1856. 5-3754 is the 
brightest and thus best qualified for detailed stud- 
ies aiming at a determination of its radius, sur- 
face gravity and gravitational redshift by means 
of X-ray spectroscopy, allowing to determine the 
equation of state of matter at supranuclear den- 
sities. In this paper, we concentrate primarily on 
this object and first summarize the main obser- 
vational facts. 

RXJ1856.5-3754 (henceforth RXJ1856) was 
discovered serendipitously in a ROSAT field by 
Walter et al. 0. Using the HST Walter & 
Matthews [I| identified the X-ray source with a 
faint blue star (V ~ 26 mag) . Its distance and 
proper motion were measured with HST as well, 
yielding (117±12)pc and 0.33 arcsec/year, re- 
spectively |5lri| . With the VLT van Kerkwijk & 
Kulkarni 7 found a faint nebula surrounding the 
point source which has a cometary-like geome- 
try with a 25" tail extending along the direction 
of motion. None of the X-ray observations re- 
vealed any variability on time scales up to ten 



years. The so far best upper limit of 1.3% (2a) 
on periodic variations in the range 10~ 3 — 50 Hz 
has been established by Burwitz et al. [Hj using 
a XMM-Newton EPIC-pn observation. Chan- 
dra LETG observations with high spectral res- 
olution show a featureless spectrum that can be 
fit by a Planckian spectrum with a temperature 
of 63 ± 3 eV Compared with the optical spec- 
trum which shows a Rayleigh-Jeans slope, the X- 
ray spectrum is reduced by a factor of ~ 6. There- 
fore the overall spectrum of the source has often 
been described by a two-temperature blackbody 
model e.g. |8l9lllll3| . 

A large number of papers have been dealing 
with the nature of this compact object and the 
proposed interpretations include almost every- 
thing which has ever been discussed in this con- 
text, ranging from "normal" neutron stars with a 
stiff or soft equation of state over neutron stars 
having a quarks core to bare strange quark stars. 
In this paper, we show that the observational data 
require a normal neutron star with a rather stiff 
equation of state. 



2. THE MAGNETIC FIELD STRENGTH 
OF RXJ1856.5 3754 

The lack of any significant spectral fea- 
tures in the LETG spectrum excludes magnetic 
fields of (1.3-7) x 10 n G (for electrons) and 
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Table 1 

Parameters of X-ray-dim isolated neutron stars. 



Source Name 


P 


P 


L x kT BB 


d 


Opt. 




(s) 


(ss- 1 ) 


(ergs -1 ) 


(eV) 


(PC) 


(mag) 


RX J0420.0-5022 Q 


3.45 




2.7xl0 30 


44 


100 


B>25.5 


RX J0720.4-3125 


8.39 


(3-6)xl0~ 14 


2.6xl0 31 


85 


100 


B = 26.6 


RX J0806.4-4123 


11.37 




5.7xl0 30 


95 


100 


B>24 


1RXS J130848.6+212708 


10.31 


<6xl0~ 12 


5.1 xlO 30 


90 


100 


m 50CCD = 28 - 6 


RX J1605.3+3249 






l.lxlO 31 


92 


100 


B>27 


RX J1836. 2+5925 






5.4xl0 30 


43 


400 


V>25.2 


RX J1856. 5-3754 






1.5xl0 31 


G3 


117 


V = 25.7 


1RXS J214303.7+065419 






l.lxlO 31 


90 


100 


R>23 


a Period and temperature from XMM-Newton data (Haberl et al. 


in preparation) . 



(0.2 - 1.3) x 10 14 G (for protons), see Burwitz et 
al. |10|. This leaves the possibility of a low mag- 
netic field characteristic for millisecond pulsars 
or a high magnetic field typical for normal pul- 
sars open. Unfortunately, due to the absence of 
a periodicity the usual estimate of the magnetic 
field of RXJ1856 based on the rotating dipole 
model is not possible. Using phenomenological 
arguments based on the very small pulsed frac- 
tion in X-rays and on a comparison with other 
objects van Kerkwijk & Kulkarni [J] have argued 
that the star has a relatively low magnetic field of 
a few 10 11 G which may be marginally consistent 
with the absence of proton cyclotron lines. But 
this is not the only possibility. We estimate the 
magnetic field making use of the spin-down lu- 
minosity dE/dt ~4 x 10 32 erg/s required for pow- 
ering the cometary-like emission nebula 7 and 
of the age of the star (t~5x 10 5 years) inferred 
from its proper motion and the distance to its 
likely birthplace in the Upper Sco OB associa- 
tion 0: assuming as usual the validity of mag- 
netic dipole braking we find a period of a ~ 1.8 sec 
and a magnetic field strength of ~1.1 x 10 13 G. 
We emphasise that these figures are very sim- 
ilar to those of the second brightest object of 
this kind, the pulsating source RX J0720. 4-3125 
(henceforth RXJ0720) whose spectral character- 
istics are also very similar to those of RXJ1856 
(see below). While the estimate of dE/dt may 
be considered as rather reliable, the age derived 
from the birthplace argument is not so certain. 



However, an age of t ~ 5 x 10 5 years (with an un- 
certainty of a factor of two) is fully consistent 
with what we know empirically about the cool- 
ing of neutron stars. We therefore conclude that 
the magnetic field of RXJ1856 is probably large, 
i.e. of the order of > 10 13 G. To support this, it 
will be important to exclude the alternative hy- 
pothesis of a millisecond pulsar |7I11| by means 
of a high time resolution observation with XMM- 
Newton which has already been approved. 

3. THE FEATURELESS X-RAY 

SPECTRUM OF RXJ1856.5 3754 

The main puzzle of RXJ1856 is the obser- 
vational fact that its X-ray spectrum is com- 
pletely featureless. It has been pointed out by 
Burwitz et al. [HUD] that nonmagnetic photo- 
spheric spectra assuming a pure iron composition 
are incompatible with the measured spectrum be- 
cause the predicted Fe-L features are not detected 
with high significance. Even a solar composition 
model with its small abundance of metals leads 
to unacceptable spectral fits. Doppler smearing 
of the spectral lines due to fast rotation does not 
wash away completely the strongest spectral fea- 
tures |12I13| . On the other hand hydrogen or he- 
lium photospheres can be excluded, because they 
over-predict the optical flux by a large factor J3j . 
Therefore any nonmagnetic photosphere can be 
firmly excluded. 

As a remedy it has been proposed that the star 
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Figure 1. Blackbody fits to the optical and X-ray spectra of RX J1856. 5-3754 for a two-component model 
(a) and a model with a continuous temperature distribution (b). 



has no atmosphere but a condensed matter sur- 
face |10I15| . Such a surface is expected to be re- 
flective in the X-ray domain |16I17) which could 
help to explain the low X-ray/optical flux ratio. 
Condensation of surface matter requires low tem- 
peratures and strong magnetic fields. To con- 
dense hydrogen at a temperature of kT = 63 eV a 
magnetic field of 5 x 10 13 G is required jTHj. For 
iron it is not clear whether a condensate can exist 
at all. According to Lai ^B] the cohesive energy 
of iron is uncertain, but condensation may oc- 
cur at 3 x 10 11 G (at kT = 63 eV) while Neuhauser 
et al. {H^ conclude that iron cannot condensate. 
Another problem is that the optical properties 
(the reflectivity depending on photon energy, po- 
larisation and magnetic field angle) of a con- 



densed matter surface have not been calculated 
taking into account the effects of atomic energy 
levels and of the solid state structure. Simpli- 
fied models which calculate the reflectivities of a 
highly magnetized, high density electron plasma 
have been performed by Lenzen & Triimper |16) . 
Brinkmann |17|. Zane et al. [201 an d Turolla et 
al. [TJ]. They yield X-ray reflectivities of typi- 
cally 10 — 30 % which smoothly depend on fre- 
quency. Close to the cyclotron frequency ~ 50 % 
is reached. We note that this is not sufficient to 
fully explain the reduced X-ray flux which would 
require a reflectivity of ~ 85 %. 

Another possibility to suppress atomic line fea- 
tures is to employ super-strong magnetic fields. 
The energy levels for iron atoms in strong mag- 
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netic fields (10 13 G) calculated by Neuhauser et 
al. |19| are pretty much smeared over the en- 
tire LETG energy range and have a spacing of 
50 — 100 eV. That could be easily resolved by the 
LETG (resolution < leV). However, the individ- 
ual levels show a B -dependence and therefore 
a smearing is expected to take place if the flux 
is integrated over the whole stellar surface. For 
a dipolar field with a factor of two variation of 
the magnetic field between pole and equator the 
lines would be broadened by 80 — 300 eV . Thus it 
is plausible, that the combination of a dense level 
structure of the magnetic atoms with a dispersion 
of the magnetic field produces a spectrum which 
appears as a continuum seen with the LETG. We 
believe that this is the most promising model for 
explaining the featureless X-ray spectrum, but it 
has to be confirmed by detailed calculations. 

4. THE ABSENCE OF PERIODIC 
VARIATIONS OF RXJ1856.5 3754 

The upper limit of 1.3% for the amplitude of 
periodic variations requires a pretty good align- 
ment between the rotational axis and the line of 
sight if the overall spectrum is represented by a 
two-temperature blackbody model |12j . We note 
that this constraint can be somewhat (but not 
fully) relaxed if the X-ray flux were reduced due 
to reflection effects because the size of the X- 
ray emitting spot would be increased. Another 
possibility is that the rotational frequency of the 
neutron star is larger than the present limit of 
observations (50 Hz). This will be checked by a 
forthcoming accepted XMM-Newton observation 
with the EPIC pn-CCD camera in using its high 
time resolution mode. 

5. COMPARISON WITH RELATED 
OBJECTS 

(1) RXJ0720.4-3125: The second brightest of 
the sources listed in Table 1 is RXJ0720 which 
looks like a twin of RXJ1856 in many respects: It 
was discovered in the ROSAT all sky survey and 
has a (blackbody) temperature of kT = 79 ± 4 eV. 
The X-ray source was optically identified with 
a B = 26.5 mag object by Motch & Haberl |2T] 



and Kulkarni & van Kerkwijk |22j. The opti- 
cal to X-ray flux ratio is about 5. RXJ0720 
shows X-ray pulsations with a period of 8.3914 sec 
and a pulsed fraction of ~10% [SHI- A timing 
analysis including all available archival and new 
ROSAT, Chandra, and XMM-Newton data indi- 
cates that the source is spinning down with a rate 
of dP/dt ~ 10~ 14 ss _1 24 . When interpreted in 
terms of magnetic braking this results in a mag- 
netic field of B-2 x 10 13 G. 

A broad spectral feature at ~ 270 eV was found 
in the EPIC-pn-spectra and when interpreted 
as a proton or iron cyclotron resonance yields 
a magnetic field strength of ~ 5 x 10 13 G and 
< 2.5 x 10 13 G, respectively, consistent with the 
value derived from the spin down |25| . 

(2) 1RXS J130848+2127 (=RBS1223): A 
strong broad band absorption feature was de- 
tected in the spectrum of this 10.31 s pulsar, with 
a similarly inferred magnetic field strength of 
(2 - 6) x 10 13 G ESJ. This pulsar shows a double 
peaked pulse profile. 

The similarities of these three sources 
strengthen the case of a strong magnetic field 
of RXJ1856. Actually the main difference be- 
tween them may be in the configuration of spin 
axis, magnetic axis and the line of sight. 

6. MODIFICATION OF THE 
PLANCKIAN SPECTRUM 

Yet another similarity between these sources 
regards the shape of the X-ray continuum. Al- 
though in the cases of RXJ1856 and RXJ0720 
a Planckian distribution yields a very good fit 
to the data a modified Planckian of the type 
E a x Planckian fits even somewhat better. This 
has been first discussed by Burwitz et al. JU| f° r 
RXJ1856 for which a = 1.28 ± 0.30, but also ap- 
plies to RXJ0720 which has a = 0.98 ± 0.30. In 
compensation for the steeper continuum the in- 
terstellar column densities njj and temperatures 
of RXJ1856 are reduced by ~ 50 and ~ 10 %, re- 
spectively. We note that such a broad band mod- 
ification does not affect the discussion of the ab- 
sence of (narrow band) line features. But it may 
be connected with an energy dependent reflection 
coefficient of the emitting surface. 
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Figure 2. The mass-radius relations for various equations of state for the nuclear matter according to 
The thick dashed curve (red) represents the apparent neutron star radius derived from both the 
two-component and continuous temperature blackbody models. 



7. A LOWER LIMIT TO THE RADIUS 
OF RXJ1856.5 3754 

In the absence of any spectral features the de- 
termination of the neutron star radius hinges on 
the distance, for which we adopt d = 117 pc |5I6| . 
In the following we first use the parameters of a 
two component blackbody model for the optical 
and X-ray spectrum of RXJ1856 (Burwitz et al. 
8 ) which is shown in Fig. la. The blackbody 
radius and temperature of the X-ray emitting 
hot spot derived from the Chandra LETG spec- 
trum are R x = 4.4km and kT x = 63 eV, re- 
spectively. The optical spectrum is interpreted 
as the sum of the Rayleigh-Jeans spectra of 
both the hot and the cool component. This fixes 
(R opt ) 2 x T opt + (R x ) 2 x T x = 33 eV x (17 km) 2 . 
The condition that the optical spectrum of the 
cool component does not show up as a deviation 
in the X-ray spectrum limits the corresponding 
temperature to kT opt < 33 eV at the 3cr level [5]. 
Using these figures we find for the radius of the 
neutron star R= (R 2 t + R|)V2 > 16.5 km (3cr). 



As an alternative we discuss a model with a 
continuous temperature distribution (c.f. Fig.^j 
of the form 



T = T 



hot 



1 



(1) 



The best fit to both the optical and the X-ray 
spectrum yields a central temperature of the hot 
spot Tj^i = 82 eV, an angular size of the hot spot 
9q = 40° and 7 = 2.1. In this case the neutron star 
radius is 16.8 km. We note that such hot spots 
may be caused by an anisotropic heat transport in 
strong magnetic fields (> 10 13 G [27|) or by polar 
cap heating in a millisecond pulsar 

The quoted apparent radii R= 16.5 km for the 
two-component model (Fig.^,) and R= 16.8 km 
for the continuous temperature distribution 
model (Fig.^>) represent rather lower limits since 
they have been derived under the assumption of 
blackbody emission. These apparent radii R are 
related to the true stellar radius Ro which is usu- 
ally quoted in the literature by 



R = R (1-Rs/R r 1/2 



(2) 
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where R s = 2GM/c 2 is the Schwarzschild ra- 
dius. For a standard neutron star of 1.4 solar 
masses the true radii are R = 14.0km (Fig.QJi) 
and Ro = 14.1km (Fig-QJ)), respectively, and thus 
considerably larger than the usual standard ra- 
dius of 10 km. This implies a rather stiff equation 
of state. We note that the same conclusion was 
reached by Braje & Romani |T2] using a two com- 
ponent model and similar arguments. In order to 
compare our results in more detail with the pre- 
dictions of theoretical neutron star models we use 
the mass-radius diagram given by Pons et al. 9 . 
This diagram is shown in Fig.|2to which we have 
added a curve corresponding to the apparent ra- 
dius derived in this paper. It is evident that the 
lower limit represented by the thick dashed curve 
excludes the quark star discussed in |2J with high 
significance. It even excludes the quark star mod- 
els discussed by Schertler et al. [2Hj which have 
even smaller radii, as well as neutron stars with 
strange quark matter cores discussed in the same 
paper. In view of the small uncertainties in the 
distance (117±12)pc of RXJ1856 and the fact 
that the blackbody leads to an underestimate of 
the emitting area we conclude that this object is 
a neutron star with a stiff equation of state and 
that the possibility for a quark star and a quark 
core star can be ruled out with high confidence. 
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